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SUMMARY: High-frequency waves (5mHz to 20mHz) have previously been sug- 
gested as a source of energy accounting partial heating of the quiet solar atmosphere. 
The dynamics of previously detected high-frequency waves is analysed here. Image 
sequences are taken using the German Vacuum Tower Telescope (VTT), Observa- 
torio del Teide, Izana, Tenerife, with a Fabry-Perot spectrometer. The data were 
speckle reduced and analysed with wavelets. Wavelet phase-difference analysis is 
performed to determine whether the waves propagate. We observe the propagation 
of waves in the frequency range lOmHz to 13mHz. We also observe propagation 
of low-frequency waves in the ranges where they are thought to be evanescent in 
regions where magnetic structures are present. 
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1. INTRODUCTION 



An ongoing debate about the heating mecha- 
nism of the solar atmosphere and the role of waves in 
it (Ghosh, 2002) present the motivation for studying 
the dynamics of high- frequency waves. This work is a 
continuation of the wave analysis, parts of which are 
already presented in previous papers (Andic, 2007a, 
2007b and Andjic, 2006) and focuses on waves with 
frequencies from ImHz to 22mHz. Special attention 
will be given to interpretation of the results for the 
frequencies above lOmHz. 

The temperature of the solar atmosphere 
varies from a minimum value in the photosphere 
(around 4 • lO^K) to a maximum typically found in 
the corona (1 • lO^K). What supplies the energy nec- 
essary to form this temperature difference is still 
under discussion. Some authors claim that high- 



frequency waves carry the necessary energy (Ulm- 
schneider, 1971a, 1971b, 2003; Stein and Leibacher, 
1974; Kalkofen, 1990, 2001; Wedemeyer-Bohm et al. 
2007) while other authors claim that low-frequency 
ones are the main source (Wang, Urlich and Coroniti, 
1995; Jefferies et al. 2006). 

Recent work by Fossum and Carlsson (2005) 
and Andic (2007a) found no observational evidence 
for flux energetic enough to promote the acoustic 
heating proposed for high-frequency waves (Ulm- 
schneider, 1971a, 1971b, 2003; Stein and Leibacher, 
1974; Kalkofen, 1990, 2001). Following these ob- 
servational results Wedemeyer-Bohm et al. (2007) 
present calculations using 3D models (Wedemeyer 
et al., 2004) and state that high-frequency acoustic 
waves do in fact have a role in the energy supply to 
the corona. 

A simple assumption is that acoustic waves 
will propagate upward, form into a shock and there- 
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fore dissipate energy. This explanation assumes that 
the magnetic field is not required for the propaga- 
tion of acoustic waves. Work by Rosental et al. 
(2002) and De Pontieu et al. (2004) have shown oth- 
erwise. Rosental et al. (2002) present 2D model 
of wave propagation in the presence of a magnetic 
field, and they conclude that the presence of mag- 
netic fields significantly complicates the waves and 
their associated dynamic. Jefferies et al. (2006) 
states that waves which were previously considered 
to be evanescent (Thio, 2006) can propagate when 
magnetic fields are present. 



2. OBSERVATIONS 



The data presented here are the same as al- 
ready presented in previous papers (Andic, 2007a, 
2007b; Andjic, 2006). 

The spectral lines Fe I 543.45nm (g^, — 0) 
and 543.29nm {g^ = 0.335) were used and 2D 
spectroscopy was performed using the German Vac- 
uum Tower telescope (VTT), Observatorio del Teide, 
Izana, Tenerife, with the Fabry-Perot spectrometer. 
The data were obtained during the mornings of 22 
and 24 June 2004, with excellent seeing conditions 
and with the use of an adaptive optics system (Berke- 
feld, Soltau and von der Luehe, 2003). The solar disk 
centre was targeted and the data were obtained in 
bursts of images. The exposure time of an individual 
image was 30ms and the cadence between successive 
images was 0.25s. The time taken to scan the line 



profile was 28.4s. This gave us a Nyquist frequency of 
17.6mHz (Grenander, 1959, various chapters). The 
field of view was 38" x 20" and the data sequence ob- 
tained on 22.06.2004 (DSl) lasted one hour, while the 
data sequence obtained on 24.06.2004 (DS2) lasted 
40 minutes. The duration of the sequences was cho- 
sen to improve reliability of analysis results (Chat- 
field, 2003, various chapters). 

No data were acquired specifically in support 
of the VTT observations, so we are unable to spec- 
ify the exact state of the photospheric magnetic field 
during the observations. Nevertheless, the MDI in- 
strument on SOHO was running in a mode where full 
disk magnetograms were obtained every 96 minutes. 
Comparison of MDI with the VTT data indicates 
that the central VTT pointing was similar to MDI, 
with an error of ±10". A G-band was used to ensure 
that the correct solar features were observed, i.e. no 
bright points in DSl and an abundance of magnetic 
structures in DS2. 

The data set DSl was taken in areas where 
no visible magnetic structuring was present. The 
MDI data revealed that the region scanned within 
the VTT field of view was in a typical quiet region, 
mostly unipolar (DSl, from 22.06.2004), with a mod- 
erately intense magnetic network in the data set DS2, 
from 24.06.2004. A filtergram of the data set DS2 
shows an abundance of g-band bright structures. 

Figs. 1 and 2 reveal the magnetograms taken 
by the MDI instrument at times corresponding to 
the time of the observations. It is obvious that both 
data sets are taken near the solar disc centre without 
significant difference in the inclination angle of the 
lines of sight. 




Fig. 1. The magnetogram of the whole solar disc from the SOHO instrument, MDI. The left panel represents 
the magnetogram taken on 22.06.2004 just before the observations commenced, while the right panel shows 
the magnetogram taken on 2^.06.2004. The tips of the arrows on both panels point at the location of used 
field of view. 
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Fig. 2. The zoomed part of the Fig.l. The left panel represent data set DSl and the right panel DS2. The 
overploted white contours are curves marking a magnetic field of 10 Gauss. The red squares represent the 
approximate location of the field of view. 



Based on these data, we conclude that the 
VTT time series data set DSl corresponds to typ- 
ical quiet Sun conditions, while DS2 has part of the 
magnetic network in the field of view (Phillips, 1992). 



DATA REDUCTION AND 
ANALYSIS METHODS 



The data were subjected to the reduc- 
tion methods already explained in previous papers 
(Andic, 2007a, 2007b and Andjic, 2006). 

As a consequence of the reduction methods 
and implemented subsequent coaligment of the re- 
constructed images the field of view and duration of 
the data sets was reduced. Consequently, data set 
DSl has a field of view 32.3" x 14.4" with a dura- 
tion of 52.54min, while data set DS2 has a field of 
view 29.3" x 7.4" with a duration of 33.13min after 
reduction. 

We estimate the formation layer of the line 
cores (Philhps, 1992; del Toro Inicsta, 2003). We 
based the estimate for the formation height of the Fe 
I 543.45nm and 543.29nm line cores on the location 
of optical depth unity at the central wavelength as 
calculated from a 3-D radiative-hydrodynamic sim- 
ulation of Asplund et al. (2000). This method is 
described in detail in the work of Shchukina and Tru- 
jillo Bueno (2001). The formation heights in NLTE 
are 258.5km and 588.7km for Fe i 543.29nm and Fe 
I 543.45nm, respectively (N. Shchukina, private cor- 
respondence) . 

The wavelet analysis used here is described in 
detail in previous papers (Andic, 2007a, 2007b and 
Andjic, 2006). The following criteria are applied to 
remove spurious oscillations: 

* The signal curve is tested against spurious de- 



tections of power that may be caused by Pois- 
son noise, assuming that it is normally dis- 
tributed and follows a distribution with two 
degrees of freedom. A confidence level of 99% 
is calculated by multiplying the power in the 
background spectrum by the values of cor- 
responding to the 99th percentile of the dis- 
tribution (Torrence and Compo, 1997; Math- 
ioudakis et al., 2003). 

The signal curve is compared with a large 
number (1500) of randomised time series with 
an identical counts distribution. By compar- 
ing the value of power found in the input sig- 
nal curve with the number of times the power 
transform of the randomised series produced 
a peak of similar power, the probability of de- 
tecting non-periodic power is calculated for 
the peak power at each time-step. This infor- 
mation is used to remove all spurious power 
from our results. (Banerjee et al., 2001) 
All oscillations of duration less than 1.5 cycles 
were excluded by comparison of the width of 
the peak in the wavelet power spectrum with 
the decorrelation time. This is undertaken to 
distinguish between a spike in the data and a 
harmonious periodic component at the equiv- 
alent Fourier frequency, thus defining the os- 
cillation lifetime at the period of each power 
maximum as the interval of time from when 
the power supersedes 95% significance until 
it subsequently dips below 95% significance 
(McAteer et al., 2004). To obtain the number 
of cycles the lifetime is divided by the period. 
All oscillations with a power of less than 15% 
of the maximum power in the time sequence 
were excluded. Wavelet analysis with the 
Morlet wavelet as the mother wavelet has a 
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tendency to also detect spurious oscillations. 
Wavelet analysis is performed for the flat- 
fields, which contain noise, and the power of 
the oscillations detected is noted and com- 
pared with those found within the dataset. 
This gives us a lower limit to the power of 
the observed oscillations. 
* All oscillations with periods of more than 
788.1 s (DSl) and 496.95 s (DS2) were ex- 
cluded, since they might be due to edge effects 
arising from the finite time span of our data. 
To study the propagation characteristics of 
the detected waves (Andjic, 2006; Andic, 2007a and 
2007b), wavelet phase coherence analysis is used 
(Bloomfield et al., 2004). When the power contained 
in the time series is studied, differences between tra- 
ditional Fourier analysis and wavelet analysis are 
small. Fourier analysis, without timing information, 
has marginally different confidence levels when com- 
pared to wavelet analysis averaged over time. 

However, the phase relations between the 
wave packets can differ in time or with a change 
in the local topology. These changes will be lost in 
Fourier analysis due to the lack of timing informa- 
tion. This is mostly because phase can take negative 
values as well, while power is always positive. There- 
fore, wavelet analysis is applied in this study. In the 
work by Bloomfield et al. (2004), a detailed descrip- 
tion of this method is given. The equations used 
for the calculation of the phase difi^erence and the 
phase coherence are stated in Table 1. in the work 
by Bloomfield et al. (2004). For pure noise this pro- 
cedure yields positive coherence. Only results with a 
coherence above 0.6 can be regarded as significant. 

Phase difference and phase coherence (Davis, 
2000, various chapters) are calculated using the in- 
tensity and velocity signals from both line cores. Ve- 
locity signals were obtained from the spectral line bi- 
sector shift (Andic, 2007a, 2007b and Andjic, 2006). 
This makes it possible to form phase differences be- 
tween velocity-velocity (V-V) signals, velocity - in- 
tensity (V-I) signals and intensity - intensity (I-I) 
signals. Signals are chosen from the same spatial 
location (i.e. same pixel in the field of view). For 
V-V and 1-1 phase difference calculations, signals 
are taken from the different spectral lines, while V-1 
phase differences are calculated from the same spec- 
tral line. This means that signals used for the V-V 
and I-l phase difference calculations are separated 
from each other by the height difference of the line 
cores in addition to the constant 13.5s gap caused by 
the scan time of the observed line profiles. 

In order to compare the results with previous 
work, phase difference spectra are formed using the 
time-averaged phase differences (Fig. 3). The phase 
difference spectra are formed so that an upward 
propagating wave leads to a positive phase difference. 
The phase differences are presented in a weighted di- 
agram, where the weighting is applied per sample 
by utilising a cross-power amplitude \fP\P-2,. Fur- 



thermore, they are binned into greyscale plots that 
are normalised to the same maximum per frequency 
bin (Lites and Chipman, 1979). Fig. 4 shows the ob- 
served V-V phase lag between the velocity fluctua- 
tions of the Fe I 543.29nm and Fe I at 543.45nm lines. 
The coherence presented with each phase spectrum 
is obtained by averaging the mean time-averaged co- 
herence (Fig. 3, panel Time av. Coherence) over all 
processed signal curves. 



4. RESULTS 



In the context of periodic phenomena phase 
angle is synonymous with phase. We calculate the 
phase using as the reference point the oscillatory sig- 
nal observed in one of the spectral lines; therefore, 
the phase difference is synonymous with the phase 
angle in this work. The term phase spectra repre- 
sents distribution of the phase angle (i.e. phase dif- 
ference) arranged in a progressive series according 
to frequency. The calculated phase spectra, in both 
data sets, shows that the observed high-frequency 
waves do propagate. Fig. 3 represents the result of the 
phase analysis using the wavelets on a pair of curves 
from data set DS2. The panels marked with LCI 
Power and LC2 Power (where LC stands for light 
curve), represent the results of the wavelet power 
transforms. Lighter shaded regions correspond to an 
increased wave power. The contours present 95% 
confldence levels, while crosshatched areas mark the 
cone of influence (COI) where edge effects can be 
important (Torrence and Compo, 1997; Bloomfield 
et al., 2006). The panel marked with Cross- wavelet 
Power represents the cross-wavelet power transform 
of both curves (in this case velocity curves), LCI 
and LC2. The panel marked with Phase Differ- 
ence presents full time series wavelet phase differ- 
ence transforms as a function of time (abscissa) and 
wave frequency (ordinate). The overploted contours 
present 10% coherence exceedance levels, while the 
crosshatched area is COI. The panel marked with 
Time av. Phase Difference represents time-averaged 
phase differences (ordinate) from the panel marked 
with Phase Difference, as a function of the wave fre- 
quency(abscissa). Midpoints correspond to the mean 
temporal phase difference. The panel marked with 
Coherence is the full time series wavelet phase co- 
herence transform. Markings are the same as for the 
'Phase Difference' panel. The 'Time av. Coherence' 
represents time-averaged phase coherence (ordinate) 
of coherence as a function of wave frequency (ab- 
scissa), where the dashed line marks mean coherence. 

The 'Phase Difference' panel shows the prop- 
agation of the high-frequency waves in different di- 
rections. Lighter shades mark upward propagation. 
A detailed explanation of how this plot is calculated 
is presented in the work by Bloomfield et al. (2004). 
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Fig. 3. Result of wavelet phase coherence analysis for the data set DS2. The wavelet analysis of both signal 
curves are shown in the panels marked with LCI and LC2 Power. The panel marked Cross-wavelet Power 
shows the cross power for the both curves. The phase difference and corresponding coherence are given in 
the panels marked with the Phase difference and Coherence. Also shown in the last two panels are the phase 
difference and coherence averaged over time. 
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Figs. 4 and 5 show that the observed direction 
of the propagation is upwards. Several phase spec- 
tra were calculated to obtain as much information as 
possible about the dynamic of the studied areas. The 
following results are presented separately for both 
data sets. 



4.1. Data set DSl 

For this dataset 46512 signal curves were com- 
pared. The top panel of Fig. 4 presents the phase 
spectrum calculated between velocity fluctuations 
(V-V) at the cores of two observed lines. The phase 
difference spectrum shows that the five-minute waves 
do not propagate. The phase difference spectrum 
displays upward propagation at higher frequencies 
up to / ~ 13mHz. The coherence is very small be- 
low the ImHz, where the disregarded frequencies are 
located, above that it varies between the values of 0.7 
to 1. The phase differences for the frequencies below 
1.26mHz can be caused by the edge effects due to the 
finite time span. 

The top panel on Fig. 5 presents the phase 
spectrum calculated between intensity fluctuations 
(I-I) at the cores of two observed lines. The co- 
herence distribution varies for this phase spectrum 
much more than for the V-V spectrum. It shows 
a decrease in the frequency range 4 — 8mHz to a 
value close to 0.6. After that, we can see the positive 
phase difference increasing. This indicates that the 
high-frequency waves observed in the intensity maps 
do travel upwards. This trend is visible even for the 
waves close to the Nyquist frequency, although in the 
frequency range after 13mHz there is sudden drop in 
phase difference, even a tendency for negative phase 
differences as well. The coherence distribution in this 
range has several peaks and drops. It is significant 
to note that in areas where the coherence drops the 
positive phase difference is not dominating the spec- 
tra. 

The black line in Fig. 6 shows average veloc- 
ity distribution calculated with the V-V phase spec- 
tra. The positive velocities mark upward propaga- 
tion, while negative velocities mark downward prop- 
agation. On average there is no propagation for the 
waves in the frequency range around 5mHz. All other 
frequencies show propagation. The upward propaga- 
tion dominates this curve, despite a strong negative 
maximum near a frequency of 15mHz. After averag- 
ing over the whole frequency range we get a propa- 
gation velocity of 26.26m/s, which demonstrates the 
domination of the upward propagation. 

The black line in Fig. 7 shows the average ve- 
locity distribution calculated with I-I phase spec- 
tra. At lower frequencies we can observe that the 
waves with frequencies around the cut-off frequency 
do propagate. The upward propagation dominates 
this curve. Averaged over the whole frequency range 
the velocity is: 19.01m/s. 




Fig. 4. The V- V phase difference and correspond- 
ing coherence distribution for calculated phase. The 
top panel represents the phase difference calculated 
for the data set DSl. The other two panels repre- 
sent the phase difference from the DS2 data set. The 
bottom panel represents the phase difference for the 
locations where the bright points appear for at least 
half of the duration of a time sequence. The mid- 
dle panel represents the rest of the field of view. On 
each panel the corresponding average coherence dis- 
tribution is plotted. The dashed red lines represent 
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coherence levels of 0. 6. 




543.45nm. On average the phase spectra show that 
oscillations appear first in velocity curves and then 
in intensity curves. This trend is apparent from the 
concentration in frequency range of 8 — 15mHz. The 
rest of the phase spectra are ambiguous. 

average velocity dis:ributIon over frequencies 



Xj*jY-^----^--*'-T^^^''ilW^ 



10 15 
frequencies [mHz] 



Fig. 6. The average velocity distribution over the 
analysed frequency range calculated with V- V phase 
difference. The black line represents distribution for 
the data set DSl; while the dashed red line and the 
dotted purple line represent the data set DS2. The 
dotted purple line represents results from the parts 
of the field of view where magnetic structures in DS2 
appear for at least half of the duration of the time se- 
ries, while dashed red line represents the results for 
the rest of the field of view. 

The top right panel of Fig. 8 presents the phase 
spectrum which is calculated between the velocity 
and intensity signal (V-I) for the spectral line Fe I 
543.29nm. The phase difference is positive at most 
frequencies, indicating that signals appear first in ve- 
locity and then intensity curves. The coherence dis- 
tribution is lowest for this phase spectra. The co- 
herence distribution tends to drop in value when the 
intensity curves are analysed. This might be caused 
by the fact that intensity curves tend to have more 
spurious oscillations than velocity curves (Carlsson, 
private communication) . 

overage velocify disfribufion over frequencies 




frequencies [mHz] 



Fig. 7. Same as in Fig. 6 but calculated with /-/ 
phase difference. 



Fig. 5. Same as Fig. 4 but the /-/ phase difference 
is presented. 



Fig. 8 shows the phase spectrum calculated be- 
tween the velocity and intensity signal (V-I). The top 
left panel represents results for the spectral line Fe I 
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4.2. Data set DS2 

For this datasot 21682 signal curves wore com- 
pared. The phase differences for the frequencies be- 
low 2mHz can be caused by the edge effects due to 
the finite time span. Those frequencies are not taken 
into consideration in this analysis. Results for this 
data set tend to follow the general trend from the 
previous data set above the frequency / ~ 7mHz. 
The middle and bottom panels in Fig.4 illustrate this 
trend. However, at lower frequencies the phase dif- 
ference is approximately zero. Moreover the coher- 
ence distribution is relatively low, and varies mucli 
more than for the previous data set. The phase dif- 
ference is calculated separately for the areas where 
the bright-points appear (bottom panel) and sepa- 
rately for other locations (middle panel). Except the 
obvious lower number of curves there is no difference 
in behaviour. 

The real difference can be noticed for the av- 
erage velocity distribution. The dotted purple line 
in Fig. 6 shows the average velocity distribution cal- 
culated with V-V phase spectra for the areas where 
the bright-points appear. In the frequency range be- 
low 6mHz the average velocities show that there is 
more propagation for areas where the bright-points 
appear. Also, in the high-frequency range, above 
12mHz, there is a tendency for greater velocities. 
The dashed red line represents the average velocity 
distribution for the rest of the field. This line shows 
that the behaviour in the rest of the field is more 
similar to the behaviour of DSl (solid black line). 
Averaged over the whole frequency range the veloc- 
ity is: 74.43m/s for areas where the bright-points 
appear and 23.88m/s for the rest of the field of view. 

An average velocity distribution calculated 
with I-I also shows different behaviour in the lower- 
frequency range for different fields of view. The view 
with the intense magnetic flux concentrations (dot- 
ted purple line, in Fig.7) shows different tendencies 
for propagation with greater velocities. The rest of 
the field of view (dashed red line) shows behaviour 
similar to DSl (solid black line). Averaged over the 
whole frequency range the velocity is: 22.52m/s for 
areas where the bright-points appear and 10.47m/s 
for the rest of field of view. 

The phase spectrum for the intensity waves 
shows similar behaviour to DSl. Where the coher- 
ence is reasonably high we can observe that the waves 
do travel upward (middle panel in Fig. 5). The over- 
all the coherence is much lower than for the previ- 
ous data set. The bottom panel represents results 
for the areas where the bright-points appear. Apart 
from the obvious smaller number of curves there is 
no difference in trend. 

The phase spectrum V-I calculated for the 
spectral line Fe I 543.45nm is presented in middle 
and bottom row of Fig. 8. Here the coherence is very 
variable and lower than for the DSl. The middle and 
bottom rows of Fig. 8 show the phase spectra for V-I 
signal. The middle panels show trends that are sim- 
ilar to the top panels (DSl). For the spectral line Fe 
I 543.45nm (left middle panel of Fig. 8) it is obvious 
that waves appears first in velocity and then in inten- 
sity curves for the frequency range 8 — ISmHz. The 



rest of the spectra are more ambiguous. The spectral 
line Fe I 543.29nm (middle left panel of Fig.8) shows 
a slightly larger concentration of the phase difference 
in a positive range for the frequencies above 5mHz, 
similar to those observed for DSl (top left panel of 
Fig.8). Both show the tendency for waves to appear 
first in velocity and then in intensity curves. 

V-I phase spectra, in the right bottom pan- 
els of Fig.8, are ambiguous due to the small num- 
ber of curves. The tendency for the positive phase 
difference is clearly noticeable only above frequency 
llmHz. The left bottom panel of Fig.8 shows be- 
haviour similar to the left middle panel. 



5. DISCUSSION AND CONCLUSIONS 

The data sets used in this work were also used 
in previous papers (Andic, 2007a, 2007b and And- 
jic, 2006). The possible problems in observing high- 
frequency waves were discussed in detail in those pa- 
pers. 

To observe the wave propagation one has to 
observe the wave at two different places in space- 
time. In this work, we observed the waves that leave 
signatures of propagation in the cores of the two dif- 
ferent spectral lines. The signals we observed do ap- 
pear first in the line core of Fe l 543.29nm and then 
in the line core of Fe i 543.45nm. This should indi- 
cate wave propagation. To interpret these observa- 
tions we used the model described in Asplund et al. 
(2000). And on the basis of this model we used the 
method described in detail in Shchukina and Tru- 
jillo Bueno (2001) to obtain the atmospheric heights 
on which observed spectral lines are formed. How- 
ever, as described in several chapters of the book by 
del Toro Iniesta (2003) the assumption that a cer- 
tain spectral line is formed at a certain atmospheric 
height is valid only as much in the used atmospheric 
model. Therefore, we have to state that the inter- 
pretation and all quoted velocity values are heavily 
dependent on the model and they are correct only as 
much as the assumptions of the model. 

An additional problem is the assumption that 
the signals observed in both line cores are separated. 
Only under that assumption can the observed phase 
difference actually indicate the wave propagation. 
Unfortunately, this assumption is also based on the 
atmospheric model and the assumed independence 
of the line formation of the different spectral lines. 
Therefore, this assumption bears similar uncertainty 
as formation heights. 

In both data sets no propagation for the low- 
frequency waves - below 5mHz - was observed. This 
is expected since most of the models (Ulmschnei- 
der, 1971a, 1971b, 2003; Stein and Leibacher, 1974; 
Kalkofen, 1990, 2001; Wedemeyer-Bohm et al. 2007) 
predict that in the lower Solar atmosphere waves 
with such frequencies will be evanescent, i.e. non- 
propagating. 

This work concentrates only on the waves 
propagating in the direction parallel to the line of 
sight. This preference is initially chosen due to pre- 
vious works on the modeling of the photospheric 
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magnetic field. Subsurface magnetic fields are con- 
tinually stretched and distorted by convective flows 

[Emonet feCattaneo 2001] resulting in the emergence 
of flux bundles in the photosphere. Newly emerged 
fields arc swept towards the boundaries of granules 
and super-granules where they interact with pre- 
existing magnetic fields. In mixed-polarity regions, 
with mean fiux density of a few Gauss, the recycling 
time scale is about 40 hours [Title & Schrijver 1998] . 
The flux tubes are less dense than their surround- 
ings and therefore buoyant. This buoyancy keeps 
the flux tubes nearly vertically oriented in the pho- 
tosphere. A flux tube with a weak magnetic field 
in thermal and hydrostatic equilibrium with its sur- 
roundings is unstable against vertical displacements 
of plasma within the tube [Parker 1978] . This implies 
that in the quiet photosphere most of the magnetic 
flux tubes are vertical (i.e. parallel to the line of 
sight). 

However, recent studies have shown complex- 
ity of the magnetic fields of the photosphere and 
chromosphere. Flux tubes, which usually channel 
the oscillations upward, have different propagation 
conditions, angle of the spread and orientation at 
the heights analysed here. [Ballegooijen & Hasan 
2003, Fig. 6] Due to this complexity of the magnetic 
fields in the photosphere and chromosphere it is not 
clear what percentage of the oscillations propagating 
upward can be determined with observations along 
the line of sight is normal to the solar surface. To 
fully analyse waves propagating along the inclined 
flux tubes the information about the magnetic field 
is necessary. However, MDI does not give us informa- 
tion on the comparable resolution with our ground- 
based data. Therefore, we analysed only the waves 
propagating parallel to the line of sight. To anal- 
yse the rest of the waves one needs high-resolution 
magnetograms. 

5.1. Data set DSl 

In this data set wave propagation was ob- 
served in the frequency range from « 7 to 13mHz. 
Although Krijger et al. (2001) claims that near this 
range there is steepening into weak shocks of acoustic 



waves on the way up, the resolution of the spectral 
line in that area makes it impossible to confirm this 
conclusion in our case. Since our Nyquist frequency 
is 17mHz, close to this frequency we have only two 

points per wave and this introduces uncertainty into 
the calculated phase differences for those frequencies. 
Also, the steepening which Krijger et al. (2001) ob- 
served starts in the range where we still observe a 
clear upward propagation in V-V maps (top panel of 
Fig. 4). There is the possibility that the drop in the 
phase difference after 12mHz and subsequent reduc- 
tion in the concentration of the phase differences is 
due to an increase in the propagation velocities of 
the waves. This matter will be resolved with more 
sensitive instruments. The propagation is less clear 
in the I-I spectra (Top panel of Fig. 5), where there 
is additional phase difference concentration around 
15mHz. 

Moreover, top panel Fig. 5 shows that around 
a frequency of 4mHz there is a very slight tendency 
towards negative phase differences just after 5mHz, 
which indicates that the observed waves are travel- 
ling downwards. This is confirmed by the average 
propagation velocity distribution that shows a nega- 
tive peak in the black curve in Fig. 7. Both findings 
indicate that in this range there is a tendency for 
downward waves. However, the V-V phase difference 
in the top panel of Fig. 4 does not show such a trend. 
This might indicate that either the sources of the 
velocity and intensity signals around this frequency 
range are not the same or there is an influence of 
NLTE effects in the core of the higher spectral line, 
Fe I 543.45nm. 

The top two panels of Fig. 8 show the phase 
spectrum calculated between the velocity and inten- 
sity signal (V-I) for the spectral lines Fe I 543.45nm 
(left top panel) and Fe I 543.29nm (right top panel). 
These panels show that the phase spectrum for the 
spectral lines tend to concentrate in a positive range 
(angle values from 40 to 100). This trend is obvious 
for the frequency ranges 8 — 14mHz in the top left 
panel of Fig. 8. The top right panel does show similar 
tendencies for the spectral line Fe I 543.29nm. Both 
panels indicate that the observed waves appear first 
in velocity and then in intensity curves. 
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Fig. 8. The V- 1 phase difference and corresponding coherence distribution for calculated phase, as in Fig.^. 
Left column represents phase difference calculated for the spectral line Fe I 543.45nm and the right one for 
the spectral line Fe I 543.29nm. The top row represents results from data set DSl, while the bottom row 
represents the results of the reduced field of view, where magnetic structures in DS2 appear for at least half 
of the duration of the time series and the middle row the rest of field of view for DS2. 



5.2. Data set DS2 

Results for this data set tend to follow the 
general trends seen in the previous data set in the 
frequency range 6-13mHz (middle panel of Fig. 4). 
The same trend is visible for intensity waves where 
the frequency range 7 — 12mHz shows a clear propa- 
gation (middle panel of Fig. 5). 

Here 21682 signal curves were compared. We 
investigate whether intense magnetic flux concentra- 



tions can cause the differences noted in phase spec- 
tra. Separate phase differences are calculated for the 
areas where the magnetic flux concentrations appear 
more often than in the rest of the field of view (at 
least for a half of duration of the time series) . Phase 
difference trends for those areas are similar to the 
rest of the field of view, as shown in the bottom 
panels of Figs. 4 and 5. The only apparent dif- 
ference arises from the smaller number of the anal- 
ysed curves. The difference between areas with a 
different magnetic flux concentration is more notice- 
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able in Figs. 7 and 6. Both diagrams show that 
the average propagation velocity distribution tends 
towards greater values for the areas where the mag- 
netic flux concentrations appear more often (dotted 
purple Unc). 

The middle and bottom rows of Fig. 8 show the 
phase spectra for the V-I signal. Both show the ten- 
dency for waves to appear first in velocity and then 
intensity curves. The middle panels show trends that 
are similar to the top panels (DSl). The bottom rows 
show clear differences be cause of the smaller number 
of analysed curves. This difii'erence is especially no- 
ticeable for the right bottom panel. Where the phase 
spectra trend towrds a positive phase difference is 
clearly noticeable only above frequency llmHz. 

Magnetic flux tubes, which usually channel 
the waves upward, have different propagation con- 
ditions, angles of the spread and orientation at the 
heights analysed here (Ballegooijen and Hasan 2003, 
Fig. 6). This might explain the higher phase veloc- 
ities registered for the waves from the area where 
magnetic structures arc appearing for at least half 
of the time series duration. Work by Dc Ponticu ct 
al. (2004) confirms that inclined fiux tubes are the 
cause of the leak of p-modes energies to higher lay- 
ers of the solar atmosphere. This might result in the 
noted positive phase velocity in the frequency ranges 
where one expects evanescent waves. Also, the work 
by Berger et al. (2004) states that the magnetic field 
in the plage region is concentrated in complex struc- 
tures, which arc not generally composed of discrete 
magnetic flux tubes. This might explain why we do 
not observe typical behaviour for evanescent waves 
(Deubner and Fleck, 1989, 1990;Fleck and Deubner, 
1989; Deubner et al., 1990) in our phase spectra for 
lower- frequency waves. 

5.3. Both data sets 

Krijgcr ct al. (2001) claim that registered 
high-frequency waves do propagate and that phase 
spectra tend to lose coherence and spread for higher 
frequencies because of the steepening into weak 
shock of acoustic waves on the way up. The be- 
haviour of the phase differences analysed here follows 
this pattern (Figs. 4 and 5). However, the frequency 
ranges of the observed waves which propagate and 
which lose coherence differ slighty to those found by 
Krijger et al. (2001). We observe propagation in the 
range 6 — 13mHz for velocity signals and 8 — ISmHz 
for intensity signals. The phase differences noted 
here are also larger than in Krijger et al. (2001). The 
larger phase differences might indicate detection of 
faster waves. On the other hand, the heights that are 
sampled in this work (258.5km and 588.7km for Fe 
I 543.29nm and Fe i 543.45nm, respectively) arc dif- 
ferent to the ones sampled by Krijger et al. (2001). 
In the work by Krijger et al. (2001) the difference 
between sampled heights is approximated as 140km 
while here it is 330km. This difference in the sam- 
pled height can also allow for the larger phase angles 
detected here. 

There are differences in the ranges and be- 
haviour of the waves detected in both intensity and 



velocity curves between this work and Krijger et al. 
(2001). Nevertheless, both phase spectra, V-V and 

I-I (Figs. 4 and 5) show that high-frequency waves 
observed in this work do propagate upwards. At fre- 
quencies around 13mHz the peak in the phase differ- 
ence is noticeable (Figs. 4 and 5). This is not observed 
in previous work, to our knowledge. This same peak 
is not so clear in the middle panel of Fig. 5. In both 
figures, (Figs. 4 and 5) is obvious that above 13mHz it 
is harder to establish propagation. Similar behaviour 
is noted in the work by Krijger et al. (2001), but 
for different frequencies. Krijger et al. (2001) claims 
that this is due to steepening of high-frequency waves 
into weak shocks. Here such a statement cannot be 
confirmed. Our results give three possibilities: first, 
steepening into shocks; second, fast waves that have 
angles above 360 degrees; third, the behaviour is a 
consequence of the low signal resolution around those 
frequencies. Also, there are several possible theoret- 
ical interpretations for the noted behaviour of high- 
frequency waves observed here (Ulmschncider, 2003; 
Kalkofen, 1990, 2001; Wedemeyer-Bohm et al., 2007; 
Rosental et al., 2002; De Ponticu et al., 2004). Pos- 
sibly the matter will be resolved in the future with 
instruments with higher temporal and spatial reso- 
lution. 

Coherence distributions for each spectrum 
were calculated for the signal curves after the ap- 
plication of all conditions that were imposed to re- 
duce spurious signals. Frequencies near the Nyquist 
frequency and frequencies below ImHz (DSl) and 
2mHz (DS2) were not cut out of the phase spectra 
and coherence distributions. Nevertheless, findings 
in those frequency ranges are not conclusive due to 
the unreliability of those ranges (see Sec. 2). 

The distributions of the average velocities for 
the intensity and velocity signals, (Figs. 6 and 7) also 
confirm that there is a dominating trend for upward 
propagation. Stating finite numbers for phase ve- 
locities has to be done with caution. Except the 
influences of the solar conditions at the formation 
heights themselves (Asplund et al., 2000) there is 
also an influence caused by post-focus instrumenta- 
tion (Bendlin et al., 1992) and effects of transmis- 
sion of the telescope (Berkefeld, Soltau and von der 
Luehe, 2003). Therefore, all those factors have to 
be taken into account when discussing and analysing 
the results of phase analysis. The phase angles show 
that there is propagation, but values of the propaga- 
tion velocities are highly dependant on the assumed 
values for the height and therefore suffer all of the 
problems connected with such assumptions. How- 
ever, the values of the observed phase velocities are 
greater in the areas where magnetic structures ap- 
pear for longer than half of the time series' duration. 
The fact that propagation velocities are greater in 
the areas with significant magnetic concentrations is 
in agreement with the work by Rosental et al. (2002). 
These authors claim that, in the presence of the mag- 
netic field concentration, the propagation velocities 
of high-frequency waves will be greater. Detection of 
the propagation in the low-frequency range, for the 
same data sample from DS2, is in agreement with 
the work by Jefferies et al. (2006) and De Pontieu 
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et al. (2004). Both works claim that the presence of 
the magnetic field allows the low-frequency waves to 
propagate. 

However, since we did not have high resolu- 
tion magnetograms for the data sets, the exact in- 
fluence of the magnetic field cannot be stated. Al- 
though there arc indications that magnetic fields can 
be responsible for such behaviour, there is a possi- 
bility that such behaviour is a consequence of the 
stochastic behaviour of the wave sources (Houdek et 
al. 1999). This question will be resolved with the 
high resolution magnetograms. 

Also, there is a significant difference in the 
diagrams of the intensity and velocity signals. This 
difference can be explained with the properties of the 
Radiative Transfer Equation (summary explanation 
with various references can be found in several chap- 
ters of the book by del Toro Iniesta 2003). In short, 
from a theoretical standpoint, velocity signals shows 
waves clearer than intensity signals. However, this 
docs not mean that results from intensity arc less re- 
liable and can be ignored. The observed behaviour of 
both signals can help in improving the future atmo- 
spheric models that will deepen our understanding 
of the solar atmosphere. 

Fig. 8 shows the phase spectrum calculated be- 
tween velocity and intensity signals (V-I). All panels, 
except the bottom right one, show a clear tendency 
for waves to appear first in velocity and then after 
in the intensity signal. The top right panel (spectral 
line Fe I 543.29nm, DSl) both middle panels (spec- 
tral lines Fe I 543.45nm (left ) and Fe I 543.29nm 
(right), DS2) and the left bottom panel (spectral line 
Fe I 543.29nm, DS2) show a tendency for positive 
phase angles. The trend for angles of w 90° for the 
low-frequency range 3 — 6mHz is stated as typical 
for evanescent waves in the work of Dcubncr et al. 
(1990). In Fig. 8 the middle left panel is the only 
one where this tendency can be seen for frequencies 
around 6mHz. Therefore we cannot certify this as 
typical behaviour for evanescent waves. 

Deubner and Fleck (1990) offers V-I phase dia- 
grams for chromospheric lines which are significantly 
different to the V-I diagrams found here, indicating 
that the waves observed here do not exhibit the be- 
haviour noted for the chromospheric waves. This dif- 
ference might be caused by the heights sampled here. 
Although Deubner and Fleck (1989) sample pho- 
tospheric and chromospheric heights, the frequency 
ranges are difficult to compare with our results. As 
stated in section 2, frequency ranges below 3mHz 
might be influenced by edge effects arising from the 
finite time span of our data. This gives us only fre- 
quencies 3 — 5mHz to compare. Unfortunately there 
is no visible agreement. Clear tendencies noticed by 
Deubner and Fleck (1989), Fleck and Deubner (1989) 
and Deubner et al. (1990) that indicate evanescent 
waves are not detectable in our data. It is possible 
that we do observe propagation of such waves, as pre- 
dicted by De Pontieu et al. (2004) and commented 
on by Jefferies et al. (2006) as the possible cause of 
chromospheric heating. Comparing om results with 
works by Deubner and Fleck (1989, 1990), Fleck and 
Deubner (1989) and Deubner et al. (1990) can lead 



to the conclusion that we observe several different 
kinds of waves simultaneously, which are expected in 
light of recently developed multidimensional models 
and theoretical explanations of the solar atmosphere 
(Rosental et al., 2002; Dc Pontieu et al., 2004; Wede- 
mcyer ct al., 2004; Wedemeycr-Bohm et al., 2007). 

In Figs. 4, 5 and 8 one can notice two clear 
streaks around frequencies below 3mHz. As stated 
in section 2 those ranges might be influenced by edge 
effects arising from the finite time span of our data, 
therefore those ranges were excluded from discussion 
and analysis here. 

The observed waves propagate upward in the 
frequency range 8-13mHz. We also noted indications 
that in the lower-frequency range, below 8mHz, there 
is an indication of the propagation in areas where the 
magnetic fiux concentrations appear. 
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UDK ... 



BHCOKO(})pCKBCHTHII TRJiaCH (SmXs 

20mX3) cy npeAJioaeHH Kao H3Bop rpHjaiba THxe 
cyHiiieBe aTMOc4)epe. HHHaMHKa tux, wpejir 
xo/ijio fleTeKTOBaHHx, lajiaca je OB^e npeciaBJi- 
jena. CeKBeHii;e cjiHKa cy /i;o5iiBeHe KopHC- 
TeibeM $a5pH-nepo cneKipoMeTap Jioii;HpaH 
y HjeMaiiKOM BaKyyMCKOM-Topait TejiecKony 
(BTT) Ha 06ccpBaTopiio acji Tch^c, Hsana, 
TeHepH(|)e. YloAsmu cy hotom pe/i;yKOBaHH cneKJie 



MCTOAOM H anajiHSHpaHH BaejicTHMa. JIa ce 
yTBp/jH /ja JIH nocMaTpaHH TajiacH nponarii- 
pajy npoBCziCHa jc anajiHsa (ftasHc pasJiHKC 6a3H- 
pana na BaejiCTHMa. 3a6HJLe3ejiH cmo npona- 
raiiH]}'- rajiaca y c()peKBeHTHOM oncery o/i;:10mX3 
TO 13mX3. TaKO^e, 3a5HJi.e3eHa je nponaraitHja 
Tajiaca y (|)peKBeHTHHM once3iiMa r^e th lajiacH 
He nponariipajy, na jiOKaiinjaMa r^e cy Shjig 
npHcyTHC MarHCTHe cipyKType. 
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